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SUMMARY 

The  shape  of  a  cable  used  to  tow  a  pitot-static  probe  from  a  helicopter  is  estimated 
for  airspeeds  up  to  120  kit.  These  estimates  show  the  probe  to  be  well  outside  the  influence 
of  the  rotor  wake.  Experimental  results  obtained  from  photographs  taken  from  a  chase 
helicopter  are  used  to  resolve  an  uncertainty  in  cable  drag  coefficient  and  comparisons  are 
made  with  the  theoretical  predictions  of  overall  shape.  The  effect  of  the  assumption  of  a 
weightless  cable  and  of  reduced  probe  drag  are  also  shown. 
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ABSTRACT 


Tlw  .shape  of  a  cable  used  lo  low  a  pilot-sialic  probe  from  a  helicopter  is  estimated 
for  airspeeds  up  lo  120  kn.  These  estimates  show  the  probe  to  be  well  outside  the  influence 
of  the  rotor  wake.  Experimental  results  obtained  from  photographs  taken  from  a  chase 
helicopter  are  used  to  resolve  an  uncertainty  in  cable  drag  coefficient  and  comparisons  are 
made  with  the  theoretical  predictions  of  overall  shape.  The  effect  of  the  assumption  of  a 
weightless  cable  and  of  reduced  probe  drag  are  also  shown. 
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NOMENCLATURE 

a  (krsi/wy- 

Cn  drag  coefficient  of  probe  (based  on  cross-sectional  area  of  main  body, 
i.e.  I  mf1)  in  direction  of  flow  along  probe  axis 

Cdc  drag  coefficient  of  cable  (based  on  cable  reference  diameter)  in  lateral 
direction  of  flow 

d  probe  reference  diameter  (i.e.  of  main  body) 

dc  cable  reference  diameter  (i.e.  of  vinyl  tubing) 

k  Jp-t  (]w/2)C» 

kc  \p.l  <lc  C/„ 

P  upper  end  of  cable  element 

Pa  upper  end  of  cable  length  (i.e.  suspension  point  at  helicopter) 

Pi  lower  end  of  cable  length  (i.e.  suspension  point  at  probe) 

s  distance  along  cable  of  P  from  Pa 

S  .si si 

T  tension  force  at  P 

u  horizontal  velocity  of  air  relative  to  cable  and  probe  (positive  in  x  direction), 

i.e.  helicopter  airspeed 

V  u  I  a 

iv  cable  weight  in  air  per  unit  length 

W  probe  weight  in  air 

at  cable  axes  with  origin  at  Pa.  c  axis  vertical  (positive  downwards),  and  .v  axis 

horizontal  in  helicopter  fuselage  plane  of  symmetry  (positive  rearward) 

.v,  r  linear  displacement  coordinates  in  cable  axes 

A ,Z  V/.M...-/M 

a  probe  sweepback  angle,  arctan  (a  i/r i) 

8  denotes  increment  in  quantity  prefixed  (see  Fig.  I) 

t>  angular  displacement  of  cable  at  P  from  c  axis  (increases  positively  towards 

a  axis) 

/i  ii’ai/M' 

i',i  kinematic’ viscosity  of  air 

P. i  air  density 

r  T If' 


Subscripts 

0.  1 


value  at  P».  Pi 


1.  INTRODUCTION 


In  order  to  validate  a  mathematical  mode!  of  the  Royal  Australian  Navy's  Sea  King  Mk.50 
helicpter  [i],  flight  trials  [2]  were  conducted  during  1979.  For  some  of  these  trials,  a  pitot-static 
probe  was  towed  from  the  helicopter  so  that  airspeed  measuring  devices  on  the  helicopter 
could  be  calibrated  [3].  The  probe  is  of  similar  design  to  that  developed  by  ARDU  (Aircraft 
Research  and  Development  Unit,  RAAF,  Salisbury,  South  Australia).  It  was  originally  developed 
by  the  US  Navy  for  the  purpose  of  giving  stable  flight  when  towed  from  a  helicopter,  but  further 
modifications  were  made  by  ARDU  to  improve  stability.  For  accurate  measurements,  it  was 
important  for  the  probe  to  be  well  outside  the  influence  of  the  rotor  wake.  To  ensure  this,  the 
shape  of  the  cable,  and  hence  position  of  the  probe  relative  to  the  helicopter,  were  estimated 
theoretically.  These  estimates,  together  with  the  two-dimensional  analysis  on  which  they  are 
based,  are  first  presented  in  this  document.  In  two-dimensional  steady  state  cable  problems 
where  gravitational  and  fluid  forces  are  considered,  analytical  solutions  are  possible  [4],  provided 
vertical  displacement  between  the  cable  ends  is  known.  Here,  because  cable  length,  rather  than 
vertical  displacement,  is  known,  the  solution  can  only  be  obtained  by  numerical  integration  of 
the  differential  equations  with  cable  length  as  the  independent  variable.  In  the  analysis,  there 
was  a  fair  degree  of  uncertainty  in  the  value  of  the  cable  drag  coefficient,  Coe ■  While  this  un¬ 
certainty  was  considered  unlikely  to  alter  the  conclusion  that  the  rotor  wake  did  not  affect 
the  probe  measurements,  it  was  felt  that  some  validation  of  the  estimates  was  desirable  and  would 
enable  more  confident  predictions  if  required  subsequently.  Hence,  measurements  were  made 
from  photographs  taken  from  another  helicopter  flying  alongside,  and  these  were  used  to  deter¬ 
mine  a  value  of  CV  which  gave  the  best  overall  agreement  between  theoretical  and  experimental 
results.  Finally,  the  effect  of  the  assumption  of  a  weightless  cable,  as  in  Reference  5,  and  of 
reduced  probe  drag  are  shown. 


2.  MATHEMATICAL  ANALYSIS 

Numerous  analysis  of  lowed  body  configurations  have  appeared  in  the  literature.  One 
such,  which  can  conveniently  be  utilised  in  the  present  application  is  given  in  References  6  and  7, 
in  which  a  three-dimensional  steady  state  mathematical  model  of  a  cable  used  to  suspend  a  sonar 
body  from  a  helicopter  is  derived.  The  model  may  be  applied  here,  where  the  type  of  solution 
required  is  two-dimensional  for  a  uniform  fluid  velocity.  In  References  6  and  7,  the  vertical 
displacement  of  the  lower  end  of  the  cable,  ri,  is  assumed  known  and  vertical  displacement, 
is  therefore  used  as  the  independent  variable.  Here,  total  cable  length,  .vt,  is  known  so  that 
displacement  along  the  cable,  w  is  defined  as  the  independent  variable.  Considering  an  element 
of  the  cable  at  R(see  Fig.  I ),  the  basic  equilibrium  equations  are  given  by  (see  Equations  (I  8a, A, r) 
of  Ref.  6) 


ic  sin  l>  k,  ii  ii  cos-  0 
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u  cos  0 
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and  the  cable  shape  is  defined  by  the  iclations 
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Dimensionless  qualities  are  now  introduced  as  in  References  6  and  7,  except  that  n  is 
replaced  by  ji,  giving  /r  =  ws\jW,  X  =  .v/.vi,  Z  =  r/si,  S  =  sjs\,  r  =  T/W,  and  V  =  uja,  where 
a  =  (kcSijW)  '.  Equations  (I)  to  (4)  may  be  written  in  terms  of  these  quantities  as 


tIO  ft  sin  0  —  u\u\ cos2 (I 

(5) 
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Unlike  the  case  where  Z  is  the  independent  variable,  an  analytical  solution  of  Equations  (5) 
and  (6)  is  not  possible.  Hence,  given  0  -  0t  and  r  -  T)  at  Pi,  Equations  (5)  and  (6)  are  solved 
using  the  Runge-Kutta  second  order  method.  The  displacements  X  and  Z  are  calculated  from 
Equations  (7)  and  (8)  using  Simpson's  formula  following  each  step  of  the  Runge-Kutta  method 
used  to  obtain  0  and  r.  As  described  in  Appendix  A  of  Reference  6,  this  procedure  is  particularly 
suited  to  the  Runge-Kutta  second  order  method  used,  since  the  mid-interval  values  required  by 
Simpson’s  formula  are  directly  provided  by  the  Runge-Kutta  method.  It  was  found  that  an 
integration  interval  in  S  of  0  02  gives  sufficient  plotting  accuracy  for  the  curves  presented  in 
the  figures. 

The  boundary  values,  ffi  and  ti,.  arc  given  on  resolving  in  the  x  and  z  directions  at  Pi 
(see  Fig.  I),  i.e. 
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3.  CABLE  AND  PROBE  PARAMETERS 

The  dimensions  and  weights  required  in  the  theoretical  analysis  for  both  the  cable  and 
probe  are  given  in  Table  I.  Diag  coefficients  weie  detivcd  theoretically  for  the  cable  and  experi¬ 
mentally  from  wind  tunnel  tests  for  the  piobe  as  follows. 

(a)  Cable  drag  coefficient 

The  cable  consists  of  a  7  7  stainless  steel  control  cable  of  diameter  2-38  mm 

(0  00781  ft)  and  two  lengths  of  vinyl  tubing  of  diameter  5  s  mm  (0  018  ft).  The  three 
lengths  are  taped  together  (see  Fig.  I)  with  the  steel  cable  fac  ng  the  airflow.  Ideally, 
the  drag  coefficient  of  the  resulting  section  is  expected  to  be  between  that  experienced 
for  an  elliptic  section  with  chord  to  thickness  ratio  of  2  0  and  that  t.  a  smooth  circular 
cylinder  For  airspeeds  in  the  lange  20  to  120  kn,  the  corresponding  Re  ”olds  number* 
range  is  from  3-9  I0:i  to  2-3  10',  which  is  subcritical.  From  Reference  8,  the  drag 

coefficients  for  the  elliptic  section  and  circular  cylinder  are  0-6  and  I  I  respect'vely. 
However,  there  is  likely  to  be  some  twisting  of  the  taped  tubes,  thus  increasing  tu; 
effective  diag  coefficient  (based  on  a  fixed  reference  diameter).  To  allow  for  this,  a  range 
in  drag  coefficient.  Ci><.  of  between  0  9  and  I  -5  is  assumed. 


I  452 


Reynolds  number  ml  i  i.  whe:e  i  i 


10  *  m-  s  (1-563  10  1  ft-/s). 


TABLE  1 


Cable  and  Probe  Dimensions  and  Weights 


Quantity  j 

Value 

Cable  reference  diameter,  dc 

5-5  mm  (0-018  ft) 

Cable  weight  per  unit  length,*  iv 

0-588  N/m  (0-0402  Ibf/ft) 

Cable  length,  si 

45-7  m  (150  ft) 

Probe  reference  diameter,  d  j 

38-1  mm  (0-125  ft) 

Probe  weight,  W  I 

64-1  N  (14-4  Ibf) 

Length  of  probe  ahead  of  cone 

1- 14  m  (3-73  ft) 

Cone  base  diameter 

0-305  m(l  -00  ft) 

Cone  vertex  angle 

90° 

Diameter  of  holes  on  cone 

19-1  mm  (0-0625  ft) 

Distance  of  centre  of  hole  from  cone  edge:f 
(i)  outer  row 

25-4  mm  (0-0833  ft) 

(ii)  inner  row 

38- 1  mm  (0-125  ft) 

*  Weight  per  unit  length  of  steel  cable  =  0-234  N/m  (0-0160  Ibf/ft) 
Weight  per  unit  length  of  vinyl  tubing  =  O'  177  N/m  (0  0121  Ibf/ft) 
w  0-234  t  2  -  0-177  =  0  588  N/m  (0  0402  Ibf/ft) 


+  Two  rows  of  equally  spaced  holes  (22  for  outer  row;  18  for  inner  row). 


(b)  Probe  drag  coefficient 

Photographs  have  been  taken  of  the  probe  when  freely  suspended  by  a  short  length 
of  cable  in  the  ARL  2-7  ■  2 •  I  m  (9  •  7  ft)  low  speed  wind  tunnel  (see  Fig.  2).  The 
deviation.  <h,  of  the  cable  from  the  vertical  at  the  probe  was  measured  from  the  photo¬ 
graphs.  At  all  airspeeds,  the  probe  was  tilted  at  about  1°  nose  down.  In  the  calculations 
below,  this  angle  was  assumed  zero.  The  probe  drag  coefficient,  Co,  may  be  calculated 
using  Equation  (I  I),  where  A-  —  \p  i  (W'-)C/>  and/>.i  =  I -226  kg/m3  (0- 002378  slug/ft3). 
Hence, 


H'  tan  (>\ 
Ip  i(W'-)w- 
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,o>  “m.*‘ 


for  u  in  knots 


Because  of  the  possibility  of  fitting  a  different  probe  tail  (sec  below),  it  was  decided  to 
base,  probe  drag  coefficients  on  the  ct oss-sectional  area  of  the  main  body,  which  would 
be  unaffected  by  any  such  change,  rather  than  on  that  of  the  tail  cone  base.  The  results 
are  given  in  Table  2.  A  mean  value  for  Co  of  44  0  is  assumed,  which  is  equivalent 
to  a  value  of  0  69  based  on  the  cone  base  area. 


TABLE  2 


Probe  Drag  Coefficient  from  Wind  Tunnel  Tests 


Photo  no. 

u  (kn) 

(deg) 

Co 

1 

41-3 

12-3 

44-2 

2 

51  -4 

18-8 

44  5 

3 

70  7 

32-6 

44  2 

4 

96-5 

50-8 

45-5 

5 

1 17- 1 

59  7 

43  I 

3 
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I 
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4.  RESULTS 

Figure  3  shows  the  predicted  cable  shapes  for  airspeeds  of  40,  80,  and  120  kn  and  for  values 
of  Cdc  in  the  range  0-9  to  I  j.  Also  shown  are  estimates  of  the  forward  boundary  of  the  rotor 
wake  at  airspeeds  of  20,  40  and  80  kn.  In  the  absence  of  suitable  experimental  data,  the  bound¬ 
aries  were  obtained  from  simulations  using  ARL’s  Sea  King  mathematical  model  [I],  In  the 
model,  the  rotor  wake  vortex  distribution  is  approximated  by  a  straight  elliptic  cylinder  as  given 
in  Reference  9..  For  the  simulations,  the  rotor  radius  is  9.46  m  (31.0  ft)  and  the  helicopter  weight 
is  assumed  to  be  84-6  kN  (19  000  ibf).  It  can  be  seen  in  Figure  3  that,  even  for  the  highest  value 
of  Cdc  assumed,  the  probe  is  well  outside  the  rotor  wake  for  airspeeds  between  20  and  120  kn. 

To  determine  a  value  of  Coe  more  precisely,  measurements  were  made  from  photographs 
taken  from  another  helicopter  flying  alongside  (see  Fig.  4).  The  cable  and  probe  were  observed 
to  remain  fairly  steady  in  flight  with  only  very  small  cable  oscillation.  When  the  photographs 
were  taken,  the  main  consideration  was  to  confirm  that  the  probe  was  well  clear  of  the  influence 
of  the  rotor  wake  and  was  steady  in  flight.  Accurate  determination  of  the  cable  shape  was  a 
secondary  consideration.  Therefore,  no  special  care  was  taken  to  minimize  parallax  errors. 
On  viewing  the  position  of  the  horizon  and  orientation  of  the  helicopter  on  the  photographs,  it 
was  considered  that,  in  general,  tilt  in  the  vertical  plane  was  likely  to  be  larger  than  that  in  the 
horizontal  plane.  Provided  the  probe  sweepback  angle  is  not  too  large  (i.e.  not  close  to  90°), 
this  will  generally  result  in  larger  parallax  errors  in  the  vertical  plane.  Hence,  error  estimates 
were  first  made  for  the  vertical  plane,  and,  since  these  showed  that  the  corrections  required 
were  mostly  insignificant,  corresponding  estimates  were  not  made  for  the  horizontal  plane.  In 
the  worst  case,  vertical  tilt  angle  was  estimated  to  be  9",  which,  for  a  vertical  angular  displace¬ 
ment  of  the  image  of  the  cable  of  I5n,  resulted  in  a  parallax  error  of  3  •5%.  Errors  arising  from 
lens  distortion  were  not  considered  significant.  Airspeed  indicated  by  the  aircraft  instrumentation 
was  corrected  from  measurements  taken  from  the  trailing  probe  (3).  Figure  5  shows  a  comparison 
between  experimental  (i.e.  from  the  photographs)  and  theoretical  values  of  probe  sweepback 
angle,  a.  The  theoretical  values  arc  given  for  values  of  Cnc  between  0-9  and  I  -5.  The  low  experi¬ 
mental  value  at  42  kn  is  possibly  due  to  an  unusually  significant  parallax  error  in  the  horizontal 
plane,  which  could  have  occurred  because,  when  the  photograph  was  taken,  the  chase  helicopter 
was  positioned  well  behind  the  helicopter  with  the  probe.  Hence,  allowing  for  this  uncertainty, 
a  value  for  Co,  of  I  -  3  was  chosen  as  giving  the  best  overall  agreement  between  theoretical  and 
experimental  results.  Using  this  value,  a  comparison  between  theoretical  and  experimental  values 
of  cable  angle  at  the  helicopter,  (f<>,  and  at  the  probe,  (fi,  is  shown  in  Figure  6,  where  agreement 
is  observed  to  be  good.  Figure  7  shows  the  corresponding  comparison  for  cable  shape. 

In  Reference  5,  estimates  of  cable  shape  arc  given  for  a  weightless  cable,  thus  enabling  an 
analytical  solution.  In  Figure  8,  it  can  be  seen  that  for  the  system  described  here  there  is  a  signi¬ 
ficant,  but  fairly  uniform,  difference  in  cable  shapes  obtained  with  and  without  this  assumption. 
Since  the  assumption  gives  results  showing  the  probe  to  be  closer  to  the  rotor  wake  boundary, 
weightless  cable  results  could  be  used  as  a  conservative  guide  to  the  likelihood  of  interference 
of  the  rotor  wake. 

Before  undertaking  the  present  study,  the  possibility  of  replacing  the  probe  tail  cone  with 
a  cruciform  finned  tail  was  considered.  The  purpose  of  this  was  to  reduce  the  probe  drag  so 
that  the  cable  would  remain  closer  to  the  vertical.  Figure  9  shows  the  effect  of  reductions  in 
probe  drag  coefficient  by  factors  of  a  half  and  quarter.  It  can  be  seen,  especially  at  higher  air¬ 
speeds,  that  the  main  effect  is  to  make  the  cable  at  the  probe  closer  to  the  vertical,  while  not 
significantly  reducing  the  probe  sweepback  angle.  A  reduction  in  cable  drag  would  have  a  much 
more  significant  effect. 


5,  CONCLUDING  REMARKS 

Theoretical  estimates  of  the  shape  of  a  cable  used  to  tow  a  pitot-static  probe  from  a  heli¬ 
copter  have  been  presented  for  airspeeds  up  to  120  kn  These  estimates  show  the  probe  to  be 
well  outside  the  influence  of  the  rotor  wake.  Measurements  made  from  photographs  taken  from 
a  chase  helicopter  are  used  to  obtain  a  more  specific  value  of  cable  drag  coefficient  within  the 
range  initially  assumed  from  theoretical  considerations  alone.  Using  this  value,  overall  agree¬ 
ment  between  experimental  Oe.  from  photograph  measurements)  and  theoretical  estimates  of 
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cable  shape  is  observed  to  be  good.  Results  obtained  assuming  a  weightless  cable  show  the 
probe  to  be  closer  to  the  rotor  wake  boundary  and  hence  may  be  used  as  a  conservative  guide 
to  the  likelihood  of  interference  of  the  rotor  wake.  Replacement  of  the  probe  tail  cone  with  a 
low  drag  cruciform  finned  tail,  while  making  the  cable  at  the  probe  closer  to  the  vertical,  does 
not  reduce  the  sweepback  angle  significantly.  To  produce  a  significant  effect,  a  cable  with  lower 
drag  would  be  needed. 
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PREDICTED  CABLE  SHAPES  AND  FORWARD  ROTOR  WAKE  BOUNDARIES 


FIG.  4  PROBE  TOWED  FROM  HELICOPER;  IMAGE  OF  CABLE  IS  ENHANCED 


Airspeed,  u  (kn) 


FIG.  5  COMPARISON  BETWEEN  THEORETICAL  AND  EXPERIMENTAL  VALUES  OF  PROBE 
SWEEPBACK  ANGLE 


FIG.  6  COMPARISON  BETWEEN  THEORETICAL  AND  EXPERIMENTAL  VALUES  OF  CABLE 
ANGLE  AT  HELICOPTER,  90,  AND  AT  PROBE, 


FIG.  9  EFFECT  OF  REDUCED  PROBE  DRAG  ON  CABLE  SHAPE 
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